SCIENTIFIC 

REPpRTS 




Ghost Imaging without Discord 



SUBJECT AREAS: 
QUANTUM MECHANICS 
QUANTUM INFORMATION 
QUANTUM METROLOGY 
QUANTUM OPTICS 

Received 
11 April 201 3 

Accepted 
1 May 201 3 

Published 
15 May 2013 



Correspondence and 
requests for materials 
should be addressed to 
J.H.S. (jhs@mit.edu) 



Jeffrey H. Shapiro, Dheera Venkatraman & Franco N. C. Wong 



Research Laboratory of Electronics, Massachusetts Institute of Technology 77 Massachusetts Avenue, Cambridge, Massachusetts 
02139, USA. 

Ragy and Adesso argue that quantum discord is involved in the formation of a pseudothermal ghost image. 
We show that quantum discord plays no role in spatial light modulator ghost imaging, i.e., ghost-image 
formation based on structured illumination realized with laser light that has undergone spatial light 
modulation by the output from a pseudorandom number generator. Our analysis thus casts doubt on the 
degree to which quantum discord is necessary for ghost imaging. 

Ghost images are formed from the outputs of two photodetectors: a high spatial -resolution detector that is 
illuminated by light which has not interacted with the object to be imaged; and a bucket detector with no 
spatial resolution that is illuminated by light which has interacted with the object. Neither detector's 
output suffices to produce an image, but an image emerges by cross correlating their outputs from a sequence of 
correlated irradiance patterns^"^. The initial demonstration of ghost imaging used the entangled signal and idler 
obtained from spontaneous parametric downconversion\ hence it was thought that ghost-image formation was a 
fundamentally quantum effect. A controversy soon arose, however, when theory^"^ and experiment^' showed 
that ghost images could be produced using the pseudothermal light obtained from passing a laser beam through 
a rotating ground-glass diffuser. One side in this debate argued that lensless pseudothermal ghost imaging 
could only be explained by nonlocal two-photon interference^ while the other argued that an equivalent 
quantitative description could be developed from semiclassical photodetection theory, in which the light is treated 
classically^^'^^. Later, it was demonstrated these two disparate interpretations could coexist^^, although the dispute 
continued thereafter^^'^^. 

The crux of the preceding argument about the physics of pseudothermal ghost imaging can be distilled to the 
following: Suppose the imaging-forming system is comprised of: (1) an optical source whose output quantum 
state has a positive Glauber P-representation; (2) linear optics, such as beam splitters, free-space propagation, and 
passage through an object transparency; and (3) conventional photodetection, specifically direct detection, 
homodyne detection, or heterodyne detection. Quantum optics has long held that the measurement statistics 
produced by such an arrangement are identical when computed from quantum photodetection theory and from 
semiclassical photodetection theory^^'^^. Recently, however, a new slant on the quantum versus classical debate for 
pseudothermal ghost imaging emerged, when Ragy and Adesso showed^° that quantum correlation— measured 
by quantum discord^^"^^ — plays a significant role in Gaussian-state ghost imaging. In particular, their calculations 
revealed that the signal-to-noise ratios of both lensed and lensless ghost imaging performed with classically- 
correlated Gaussian state light— the standard source model for pseudothermal ghost imaging— shared paramet- 
ric dependence that mimicked the corresponding behavior of the total correlations, i.e., the classical correlations 
plus the quantum discord, even though this source is unentangled and, indeed, has a positive Glauber 
P-representation. 

The purpose of the present paper is not to refute any of the calculations reported by Ragy and Adesso, but 
instead to present a rather different view of classical- state ghost imaging, specifically a Gaussian pure-state 
treatment in which quantum discord cannot enter. Consequently, although one may rightly argue about the 
fundamental importance of quantum discord in Gaussian mixed-state ghost imaging, no such role can possibly 
exist for all forms of ghost imaging. 

Results 

Consider the spatial light modulator (SLM) form of lensless ghost imaging shown in Fig. and implemented by 
Venkatraman et al in Ref 26. Here a pseudorandom number generator (PRNG) applies a sequence of phase maps 
to the K X K pixels of an SLM, one map for each of M normally- incident plane-wave laser pulses that illuminate 
the SLM. The modulated light is divided into signal and reference beams by a 50-50 beam splitter, with one 
propagating over an L-m-long free-space path to a charged-coupled device (CCD) detector and the second 
propagating over an L-m-long free-space path to an object transparency. A bucket detector then collects the light 
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Figure 1 | Setup for spatial light modulator (SLM) lensless ghost imaging using a pseudorandom number generator (PRNG), a charge-coupled device 
(CCD) camera, and a bucket detector. 



transmitted through that transparency, and the ghost image is formed 
by cross correlating the sequence of bucket detector outputs with 
those from the CCD. In Ref 25, the entire collection of pseudoran- 
dom phases — M phase maps each comprised of iC X iC pixel phases — 
were taken to be independent, identically distributed, random 
variables that were uniformly distributed on [0, and an express- 
ion was presented for the co variance function of the zero -mean field 
that results after far-field diffraction. Then, invoking the central limit 
theorem to justify the assumption that this far-field illumination is a 
Gaussian random process, SLM ghost imaging was argued to be 
equivalent to pseudothermal ghost imaging. The analysis in Ref 25 
employed semiclassical photodetection theory, because, as noted 
earlier, it relies on laser light that has undergone linear- optics trans- 
formations — at the SLM and in far-field diffraction — followed by 
conventional direct detection. It is essential to the present purpose, 
however, to recast that analysis in a fully quantum-mechanical form, 
and, moreover, to pay special attention to the pseudorandom num- 
bers used to drive the SLM through its sequence of phase maps. 

In lieu of the classical SLM input and output fields consider ed in 
Ref 25, let EyJ^p.t) be the positive-frequency, photons / m^ s-units 
field operator that illuminates the SLM in Fig. 1 at transverse position 
p = (x,y) and time t, and let Eout{p,t) be the field operator that results 
at the SLM's output. The former will be taken to be in the coherent 
state |£in(p,0)j with eigenfunction 

M-l 



Ein{p,t)-- 




-mTr)e' 



-ioJot 



where: o^o is the laser's center frequency; p(t) is a baseband pulse that 
is time-limited to 0 < t < and normalized to satisfy \dt\p{t)\^ = 1; 
Tr > Tp is the pulse-repetition interval; and Nt is the state's total 
average photon number. We will assume that the SLM has square 
pixels of size dX 100% fill factor, and overall dimensions D X D, 
with the ijth pixel occupying the region Aij in the z = 0 plane, where 
the SLM is located. Assuming a lossless SLM, we then have that 

EontiP^t) =Em{p,ty , for peAij and 1 < < K, where (l)ij{t) = 



(j) ijm ^or mTr ^ t ^ mTr + Tp.O^m^M— 1 and zero otherwise, 
with { (^y^: 1 < i,j < K] being the mth pseudorandom phase map. So, 
conditioned on knowledge of the pixel phases, £out (p^t) will be in the 
coherent state |£out(Pj 0) with eigenfunction 
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for p E Aij and 1 < /, j < K. The quantum form of the Huygens- 
Fresnel principle^^'^^— together with the assumptions that the 
are independent, identically- distributed random variables which are 
uniformly distributed on [0, In] and iCis an odd integer— then repro- 
duces the semiclassical mean and covariance from Ref 25 for far- field 
diffraction of the quantum signal and reference field operators, 

^out 
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that emerge from the 50-50 beam splitter in Fig. 1, where Ey{p,t) is 
the vacuum-state field operator entering the beam splitter's free input 
port. Thus, coherent-state illumination of the SLM plus the statistics 
just assumed for the pixel phases leads to pseudothermal, i.e., mixed 
Gaussian-state, behavior for the joint state of the signal and reference 
field operators impinging on the transmission mask and the CCD in 
Fig. 1. Hence the discord analysis from Ref 20 is directly applicable. 
But a pseudorandom number generator does not produce true ran- 
domness, and this has a very significant consequence with respect to 
the role of discord in SLM ghost imaging. 

Pseudorandom numbers are produced by supplying a determin- 
istic seed number to a deterministic algorithm. Consequently, it is 
entirely fair for us to have complete knowledge of the {(l)ijm} entering 
into the preceding description of SLM-based ghost imaging. It imme- 
diately follows that the joint state of the signal and reference beams 
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Figure 2 | Spatial light modulator ghost images of the MIT logo versus M, the number of illumination patterns. 
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M = 300 M = 600 M = 1200 M = 2400 M = 4800 M = 9600 
Figure 3 | Computational ghost images of the MIT logo versus M, the number of illumination patterns. 



that illuminate the object transparency and the CCD detector is the 
coherent- state product, \Es^^^ {p\t))\ER^^^ with identical eigen- 

functions given by 

£ou.(p,f-L/c)e*l'''-''lV^^ 



Ek„ip',t)= dp 



V2 



ill 
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for / = S,R, where k = ojqIc and k = IucIcdq are the wave number and 
wavelength associated with the center frequency coq, and quasimo- 
nochromatic, paraxial propagation over an L-m-long free-space path 
has been presumed. Thus the SLM ghost image, found by cross 
correlating the outputs from the CCD and the bucket detector, arises 
from Gaussian pure-state illumination of the transparency and the 
CCD, a situation in which there is no quantum discord because 
quantum discord is a measure of entanglement in the pure-state case. 
In short, quantum discord is not essential to forming a ghost image. 
Indeed, in SLM lensless ghost imaging it plays no role whatsoever. 

Discussion 

In Ref. 15 we asserted that "a ghost-imaging configuration is 'classical' 
if its measurement statistics are correctly described by the use of 
semiclassical photodetection theory, and that the configuration is 
'quantum' if a correct description of its measurement statistics requires 



the use of quantum photodetection theory." Thus Refs. 13, 14 con- 
cluded that pseudothermal, i.e., Gaussian mixed-state, ghost imaging 
was classical, because its statistical descriptions in the semiclassical and 
quantum theories of photodetection were identical. Ragy and Adesso's 
discord analysis^°, however, suggests that there is a fundamentally 
quantum aspect to pseudothermal ghost imaging, whereas we have 
argued in the preceding section that such is certainly not the case for 
SLM ghost imaging. Our argument relied on the intrinsically deter- 
ministic nature of the pseudorandom phases applied to the spatial 
light modulator in the Fig. 1 setup, as opposed to the random nature 
of laser light propagation through a rotating ground-glass diffuser. But 
how random is laser light propagation through that difiuser? Cannot it 
be characterized by a highly complicated but nonetheless deterministic 
spatially-varying refractive index? If such is the case, then conven- 
tional pseudothermal ghost imaging's statistical characterization of 
its far- field irradiance patterns merely represents our inability to spe- 
cify that refractive-index distribution and to do the multiple-scattering 
analysis of light propagation through it. Carried to its logical conclu- 
sion, this train of thought would conclude that pseudothermal 
ghost-image formation, like SLM ghost imaging, relies on product 
coherent-state illumination and hence is discord free. We will not 
argue that point, largely because it entails the question of whether true 
randomness exists in classical electromagnetism. Instead, we devote 
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Figure 4 | Setup used to collect the ghost images in (a) Fig. 2 and (b) Fig. 3. PBS: polarizing beam splitter. HWP: half-wave plate. SLM: spatial light 
modulator. HR: high reflector. 
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the rest of this section to underscoring the deterministic nature of SLM 
ghost imaging. 

Figure 2 shows SLM ghost images of the MIT logo formed from a 
setup Hke that shown in Fig. 1 with M values ranging from 300 to 
9600; see Methods for the details. That this is truly a deterministic 
setup can be demonstrated as follows. Because the SLM in Fig. 1 is 
driven by a deterministic sequence of phase maps, there is no need to 
propagate a reference beam: free-space diffraction of the eigenfunc- 
tion Er.^ {p,t) can be computed by evaluating the diffraction integral 
in Eq. (5) and cross correlating the calculated with the 

output from the bucket detector. What results is a computational 
ghost imager, first described in Ref 25 and subsequently demon- 
strated in Ref. 27. In particular, for the MIT logo whose SLM ghost 
images are shown in Fig. 2, the setup from Ref 26 yielded the cor- 
responding computational ghost images shown in Fig. 3. Here there 
is absolutely no possibility for quantum discord to play any role, 
because there is only one optical beam used to form the ghost image. 
Note that prior to the appearance of Ref 25, computational ghost 
imaging was already known in the signal processing community as 
the structured illumination version of the single-pixel camera^^. 

In conclusion, the importance of quantum discord in quantum 
communication and computation is being widely studied, and Ragy 
and Adesso^° have shown how it may be used to understand pseu- 
dothermal ghost imaging. We, however, have shown that discord is 
not universally at work in ghost imaging. 

Methods 

Figures 4(a) and (b) show the experimental setups used to obtain the SLM-based 
ghost images shown in Figs. 2 and 3, respectively. For the SLM ghost images of Fig. 2, 
a 10-mW, 795-nm laser source was divided by a 50-50 beam splitter into signal and 
reference beams, each of which was focused onto a Boulder Nonlinear Systems SLM 
unit with 512 X 512 pixels, each 15X15 fim in size with individually addressable 
phases. We focused each beam at the SLM plane with a beam waist Wq ~ 200 fim. A 
computer was used to generate a 128 X 128 array of pseudorandom phases that we 
applied synchronously to the central 128 X 128 pixels of SLMl and SLM2 at a frame 
rate of 2.5 Hz, giving rise to phase-insensitive cross correlations. The imaging object 
was a binary transmission mask bearing the MIT logo printed onto a transparency 
slide that we placed in the signal arm at a distance L = 80 cm from SLMl. A 10- cm 
focal length lens was used to focus all the transmitted light onto a single-pixel (bucket) 
silicon detector (Thorlabs PDA55). To obtain the SLM ghost images of Fig. 2, we 
placed in the reference arm a shutterless Easier Pilot CCD camera with 1600 X 1200 
pixels, each measuring 7.4 X 7.4 fim with 12-bit dynamic range. The CCD was placed 
at the same distance L from SLM2, as shown in Fig. 4(a). The central 1200 X 1200 
portion of the CCD pixels were used for imaging and re-binned into 300 X 300 
superpixels, each consisting of 4 X 4 pixels. By computing the covariances between 
each of the CCD superpixels and the bucket detector output over M trials, we obtained 
the SLM ghost images in Fig. 2 for different M. We should note that the setup of 
Fig. 4(a) was also used to demonstrate phase- sensitive classical ghost imaging in 
Ref 26 by programming SLM2 with the exact SLMl -complementary pattern to 
generate phase- sensitive cross correlations between the signal and reference beams. 

To obtain the computational ghost images of Fig. 3 we used the setup of Fig. 4(b) 
whose signal arm is identical to the setup in Fig. 4(a) but without the reference beam 
and all the optical components in the reference arm, including SLM2. Instead, we 
used Gaussian beam propagation theory in the far field^^ to compute the fields 
deterministically from our knowledge of the computer-generated pseudorandom 
phase patterns that were applied to SLMl. By computing the covariances between the 
single-pixel detector output and the computed far-field irradiance patterns over M 
trials we obtained the computational ghost images in Fig. 3 for different M using only 
a single bucket detector. 
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